Adsorption Kinetics of a Single Polymer on a Solid Plane 



S. Bhattacharya 1 , A. Milchev 1 * 2 , V.G. Rostiashvili 1 , A.Y. Grosberg 1 ' 3 and T.A. Vilgis 1 
Max Planck Institute for Polymer Research 10 Ackermannweg, 55128 Mainz, Germany 
2 Institute for Physical Chemistry, Bulgarian Academy of Science, 1113 Sofia, Bulgaria 
3 Department of Physics, University of Minnesota, Minneapolis, MN 55455, USA 

We study analytically and by means of an off-lattice bead-spring dynamic Monte Carlo simulation 
model the adsorption kinetics of a single macromolecule on a structureless flat substrate in the 
regime of strong physisorption. The underlying notion of a "stem-flower" polymer conformation, 
and the related mechanism of "zipping" during the adsorption process are shown to lead to a Fokker- 
qq | Planck equation with reflecting boundary conditions for the time-dependent probability distribution 

function (PDF) of the number of adsorbed monomers. The theoretical treatment predicts that 
■ the mean fraction of adsorbed segments grows with time as a power law with a power of (1 + 

£\J ' u ) where v ~ 3/5 is the Flory exponent. The instantaneous distribution of train lengths is 

predicted to follow an exponential relationship. The corresponding PDFs for loops and tails are also 
derived. The complete solution for the time-dependent PDF of the number of adsorbed monomers 
is obtained numerically from the set of discrete coupled differential equations and shown to be in 
perfect agreement with the Monte Carlo simulation results. In addition to homopolymer adsorption, 
we study also regular multiblock copolymers and random copolymers, and demonstrate that their 
adsorption kinetics may be considered within the same theoretical model. 

PACS numbers: 05.50.+q, 68.43.Mn, 64.60.Ak, 82.70.-y 
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I. INTRODUCTION 



CPS 

o 

The adsorption of polymers at equilibrium is fairly well understood from theoretical [l|, 0, 0, 0, 0, 0, 0] i computer 
] simulation 0, 0, 0, Op and experimental pjj [l2j points of view. On the other hand, a great deal of work exist on 
(— i . the polymer non-equilibrium adsorbtion of single polymer chains. One of the earliest MC - simulation along this line 
was implemented for single chains on the cubic lattice [HI . The authors have tested a totally irreversible adsorption 
O . model and a reversible one where a move resulting in the desorption of a segment was assigned a relative weight 
exp(xs) with Xs being the segmental adsorption energy in units fc^T. It was found that at \s > ^ksT the fraction 
, of segments in loops and trains start to deviate from their equilibrium values, i.e. the process become irreversible. 
^ • Eventually at %s ~ lOfceT the fraction of segments in loops and tails merges the corresponding values for the totally 
OO \ irreversible model. 

One of the important questions concerns the scaling of the adsorption time T a d s with the length of the polymer 
chain N. Shaffer [l4| has studied this problem using Monte Carlo simulatitions with the bond fluctuation model 
(BFM) for strong sticking energy of lOfcsT. He found that the deviation in the instantaneous fraction of adsorbed 
monomers from its equilibrium value can be described by a simple exponential decay. During the late stages, however, 
the relaxation function begins to deviate from an exponential, and the relaxation slows down considerably. According 
to Shaffer, this might be due to artefacts of the lattice model (cf.also ref.[l5j]). The main result of Shaffer [l4[ is that 
Tads ~ A 1 - 58 (for comparatively short chains A < 80). 

The same scaling has been found by Ponomarev et aZ.[l6[ who also used the BFM for N < 100. Except the energy 
gain (per segment) of e s , an activation barrier tb for a segment to desorb was introduced in this simulation, defining 
thus a "temperature" Tj, = fcgT/eb. This sets a characteristic time for the passage of a segment across the barrier 
T~b = To cxp(l/Tt,). Different adsorption dynamics has then been found, depending on the ratio of Tb and the Rouse 
time: Tb/rn ~ A -2 ^ 1 exp(l/Tf,). The case Tb/r^ <C 1 (at T& > 1 ) corresponds to strong physisorption. On the 
other hand if the chain is relatively short and the barrier is high enough (Tb is low enough) then Tb/rp/ 3> 1, which 
corresponds to chemisorption. They argue that at Tb/rji <C 1 the adsorption follows a zipping mechanism whereby 
the chain adsorbs predominantly by means of sequential, consecutive attachment of monomers, a process that quickly 
erases existing loops. In this case r a d s ~ A 157 for a SAW chain in agreement with Shaffer's results [3]. In the 
opposite limit (chemisorption), the presence of a barrier enhances loop formation in the course of adsorption. It 
was shown that even a modest local barrier discourages the tendency for zipping and switches on a new mechanism 
involving loop formation. The scaling law in that regime reads r a( js ~ A", where the exponent a = 0.8 ± 0.2. 

The irreversible chemisorption from the dilute polymer solution has been theoretically studied flTl fl8j by making 
use of the master equation (ME) method [l9( for the loops distribution function. The authors argue that the process 
is dominated by accelerated zipping when the sequential adsorption is disrupted by large loops formation. 

For strong physisorption the simple zipping mechanism, as op pos ed to the accelerated zipping, has been also recently 
considered by Descas, Sommer and Blumen [2(1. The authors [2(| used the BFM and suggested a simple theoretical 



2 



description of the corresponding adsorption dynamics based on what they call a "stem - corona" model. This leads 
to the scaling prediction, r ac i s ~ N 1+u , which is in reasonably good agreement with the simulation result. 

In the present paper we study the case of strong physisorption by means of an off-lattice dynamic MC method 
and theoretically by employing the ME-formalism. This makes it possible to describe the adsorption dynamics not 
only in terms of the average fraction of adsorbed segments but also to include train and tail distribution functions 
which furnish the main constituents of the dynamic adsorption theory. Section [TT] starts with the description of the 
adsorption dynamics model which shares many common features with the one suggested by Descas et al. (20| . Then 
we use this model within the ME-formalism to treat the time evolution of the distribution of adsorbed monomers (as 
well as the distributions of the monomers forming trains and tails) . It is shown that the problem can be mapped onto 
a drift-diffusion process governed by a Fokker-Planck equation. We obtain the numerical solution of this equation 
and discuss its consequences. In Section IIIII we briefly introduce the MC model. The MC-simulation results for 
homopolymers as well as for block- and random - copolymers in the regime of strong adsorption are given in Section 
IIVI We show that our MC-findings are in good agreement with the theoretical predictions. We summarize our results 
and conclusions in SectionfV] Some details of the train distribution function calculation are relegated to the Appendix. 



II. ADSORPTION DYNAMICS IN TERMS OF TRAIN AND TAIL DISTRIBUTIONS 

Consider a single polymer molecule (grafted with one end to a flat structureless surface) in an adsorption experiment 
which is repeated over and over again. The monomer - surface interaction is considered attractive with a sticking 
energy e = E\ — E2, where E\ and E2 are respectively the monomer energies before and after the adsorption event. 



A. Stem-flower scenario: A macroscopic law 



As indicated by earlier MC-simulation results [la, |20| , in the strong physisorption regime the process is assumed 
to follow a simple zipping mechanism. Figure QJi gives snapshots of the chain conformation as it is evident from our 
simulation results. One can see that the chain conformation can be considered within the framework of a "stem-flower" 
picture which was discussed first by Brochard-Wyart [2lj as characteristic for a polymer chain under strong stationary 
flow. Recently the "stem- flower" picture was employed in the case of non-stationary pulled polymer chain [22j. It 
should be pointed out that this picture shares many common features with the "stem-corona" model, suggested by 
Descas et al [20| . Here we reconsider it in a more systematic way and employ it as a basic model to include fluctuations 
within the ME - formalism. 

Fig. QJ> presents schematically the stem-flower scenario of the adsorption dynamics. The number of adsorbed 
monomers at time t is denoted by n(t). The nonadsorbed fraction of the chain is subdivided into two parts: a 
stretched part ("stem") of length m(t), and a remaining part ("flower") which is yet not affected by the tensile force 
of the substrate. The tensile force propagation front is at distance R(t) from the surface. The rate of adsorbtion is 
denoted as v(t) = a dn }f > , where a is the chain (Kuhn) segment length. 

A single adsorption event occurs with energy gain e and entropy loss hx{fj,a/fi2), where /Z3 and \X2 are the connectivity 
constants in three and two dimensions, respectively [23| . In result, the driving force for adsorption can be expressed 
as 

_ e- k B T\n(n 3 /n 2 ) _ F 

/drive V J 

a a 

where F = e — ksT 111(^,3/ ^2) is the change in free energy. The friction force is related to the pulling of the stem at 
rate v(t), i.e. 

fin c = Co am(t)^- (2.2) 
at 

where Co is the Stokes friction coefficient of a single bead. The equation of motion follows from the balance of driving, 
/drive, and drag force, /fric, which yields 

Co mit) = ^ (2.3) 

dt a z 

One may express m(t) in terms of n(t), if one assumes that at time t the "flower" (which is placed on average at a 
distance R(t) from the surface) is not affected by the tensile force. This means that R(t) is the size which the chain 
portion n(t) + m(t) occupied before the adsorption has started, i.e., 

a[n{t) + m(f)] v = R{t) (2.4) 
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FIG. 1: (a) Snapshots of an TV = 256 chain conformation, taken at successive time moments during the adsorption process. 
The 2-coordinate of the i-th monomer is plotted against monomer index i. (b) Stem-flower picture of the adsorption dynamics. 
The total number of adsorbed monomers at time t is denoted by n(t). The tail which, contains all nonadsorbed monomers, 
consists of a stretched part, a "stem", of length m(t), and of a nonperturbed part which is referred to as "flower". The rate of 
adsorption is v(t). The distance between the surface and the front of the tension propagation is R(t). 



where v is the Flory exponent (e.g., v = 3/5 in d = 3-dimensions) [231 ]. On the other hand, as shown in Fig. [TJ), 

am(t) ps R(t) (2.5) 

up to a geometrical factor of order unity. Therefore the relation between m(t) and n(t) is given as 

n(t) ~m{t) 1/v -m(t) (2.6) 

During most of the adsorption process the stem is sufficiently long, m(t) 3> 1, so that m(£) 5 / 3 3> m(t), i.e., 
m(t) ~ n(t) u and Eq. (|2.3p becomes 



NI/ dn(t) _ F 
a 5 



Con(<r ^T = ^ (2J) 



The solution of Eq. (|2?7|) reads 



n(t) tx 



F 

t 



-i i/(i+^) 



a 2 Co 



(2.8) 



In result, (for d — 3 where v = 3/5) one obtains a law for the adsorption kinetics, n(t) oc t - 62 , which is in a good 
agreement with MC- findings [3, 0,!!^]. In the course of adsorption the "stem" grows and the "flower" moves farther 
away from the surface. This, as it was mentioned in ref . (2p| . makes the nucleation of a new adsorption site on the 
surface less probable. 

In the late stages of adsorption the "flower" has been largely consumed and vanishes so that the non-adsorbed part 
of the macromolecule exists as a " stem" only. From this moment on the closure relation reads 

n(t) + m(t)=N (2.9) 

Comparison of Eq. (|2.9p with Eq. (|2.6p shows that this pure "stem" regime starts at n(t) > N — /V" ps N, i.e., it 
could be basically neglected for sufficiently long chains. 

The stem- flower scenario which we used in this section as well as the macroscopic equation of motion, ea. (|2.7p . are 
employed below as a starting point for the treatment of fluctuations. 



B. Time evolution of the distribution of adsorbed monomers 



Next we focus on the instantaneous number of adsorbed monomers (i.e., the total train length) distribution function 
P(n, t). The number of adsorbed monomers n and the number monomers in the nonadsorbed chain tail I are mutually 
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complementary, if one neglects the loops (we will argue below that in the strong adsorption regime the loop contribution 
is rather small and reduces mainly to loops of size unity). With this assumption, the corresponding tail distribution 
function T(l,t) reads 

T(l,t) = P(N - l,t). (2.10) 

Both P(n, t) and T(l, t) can be obtained either from the simulation or by solving a set of coupled kinetic equations. 
For the latter we use the method of the Master Equation [19( . We treat the adsorption as a sequence of elementary 
events, describing the (un)zipping dynamics while keeping in mind that within an elementary time interval only one 
monomer may change its state of sorption. Thus one can treat the (un)zipping dynamics as an one-step process, 
shown schematically in Fig. [2^i. 

In order to specify the rate constants, we use the detailed balance condition [l9( which in our case (cf. Fig. [2^) 
reads 

= e^ (2.11) 

w (n) 

where again F = e — kBT\n{^/ ^,2) is the free energy win upon a monomer adsorption event and the energy gain 
e = E\ — i?2- Detailed balance condition Eq. (|2.1ip is, of course, an approximation for the non-equilibrium adsorption 
process in question. This implies that, despite the global non-equilibrium, close to a "touch-down" point the monomers 
are in local equilibrium with respect to adsorption-desorption events. This also means that the monomer size is small 
enough as compared to the "stem" length, so that this approximation is a good one, compatible with the "stem-flower" 
picture of adsorption dynamics. 




FIG. 2: (a) Creation - annihilation of an adsorption state with n-monomers due to a single-step process. The arrows indicate 
possible single-step transitions with w + (n) and w~ (n) being the rate constants of adsorption and desorption events, respectively, 
(b) The adsorbed monomers form trains, divided by defects (loops of length unity). The total number of adsorbed monomers 
at time t is denoted by n(t). The train length, h, itself is a random number, subject to an exponential distribution D(h,t) - 
Ea. f2~32jl . 



The detailed balance requirement fixes only the ratio of the rate constants and does not fully determine their values 
which could be chosen as 

w~(n) = q[m(n)]e~ F/kBT 
w + (n - 1) = q[m{n)}. (2.12) 

In Eq. 12.121 the transmission factor q[m{n)\ is determined by the friction coefficient £ which, within our stem-flower 
model, is defined as £ = Co™- Therefore, one obtains 

q[m(n)\ = -j- = — . (2.13) 

a% a Co m 

The notation q[m(n)] implies that the stem length m depends on the total train length n and, furthermore, the 
relationship m(n) is given by the closure Eq. (|2.6|) which also holds for the instantaneous values, i.e., 

n ~m xlv -m (2.14) 
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With the rate constants from Eq. 12.121 at hand, the one-step master equation reads 

— P(n,t) = w-(n+l)P(n+l,t) + w + (n-l)P(n-l,t) 
at 

- w+(n) P(n,t) -w~{n) P(n,t) (2.15) 

or, in a more compact form 

— P(n,t) = A [w-(n) P(n,t)] + A" 1 [w + (n) P(n,t)] (2.16) 

where the finite-difference operators A, A^ 1 are defined as 

A/(n) = /(» + l)-/( n ) 
A-VW = /(n-l)-/(n) (2.17) 

The total number of the adsorbed monomers varies between 1 and N, i.e., 1 < n < N. For n = 1 the Eq. (|2.15[) 
has to be replaced by 

±P(l,t)=w-(2)P(2,t)-w+(l)P(l,t) (2.18) 

Similarly, for n = N the ME reads 

^-P(N,t) = w+(N - l)P(N - l,t) - w~(N)P(N,i) (2.19) 

Finally, the set of master equations (|2.15j) . (I2.18|) and (|2.19j) should be supplemented by the initial condition 

P(n,t = 0) = S(n-l) (2.20) 

because the adsorbtion starts from the state of a one chain end grafted at the surface. 

The equation of motion for the first statistical moment, (n) = J^n=i n -P( n iO? can De obtained from Eq. (|2.16|) by 
performing the summation by parts: 

N-l N 

9(n)A f(n) = g(N)f(N) - 5 (0)/(0) + ^ /(^)A _1 p(n) (2.21) 

n— n— 1 

where /(n) and c/(n) are arbitrary functions. Taking this into account and keeping in mind that P(N, t) = P(0, t) = 
for simplicity, the equation of motion for (n) then yields 

j t (n) = -(w-(n)) + (w + (n)) (2.22) 
With the relations for the rate constants, Eqs. (|2.12p and 12.131 this equation of motion becomes 

t m{t)±n{t) = ^- [l-e-^j (2.23) 
at a l J 

where for brevity we use the notations n(t) = (n) and m(t) — (to). The result, Eq. (|2.23[) . should be compared with 
Eq. (|2.3p derived earlier by means of a simplified physical consideration (see also paper [2(| where this result was 
obtained before us). Formally, Eq. (|2.23j) transforms back into Eq. (|2.3|) when adsorption is very weak, F/kBT <C 1. 
Importantly, Eq. (|2.23p has the same structure as Eq. (|2.3p even when the adsorption is not weak and the quantity 
F/ksT is not small; the only difference between these equations is that the effective force in Eq. (|2 . 23[) has the form 
(k B T/a) [l - e~ F l kBT \ instead of just F/a. This can be understood by the analogy with the second virial coefficient 
of interaction between the monomer and the surface. Indeed, we know that the contribution to the free energy of an 
imperfect gas due to pair collisions is proportional to the second virial coefficient rather than just interaction energy; 
similarly in the case of adsorption, the effective second virial coefficient is the quantity that describes the effect of 
monomer attraction to the wall. Thus, the zipping as a strongly non-equilibrium process can not be treated quasi- 
statically by making use of a simple "force-ballance" . The inclusion of fluctuations by employing the ME-formalism 
is important in order to obtain the correct result for the driving force. 
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1. Fokker - Planck equation and boundary conditions 

It is instructive to change now from the discrete representation, Eqs. (|2 . 16[> . (|2 . 18[) and (|2.19[) , to a continuous one, 
namely, to the Fokker-Planck equation for the distribution function Pin, t) with proper boundary conditions. This 
can be done by the substitution 

d_ 1 d 2 

dn 2 dn 2 

A-^-A + I^. (2.24) 
dn 2dn 2 { ' 

After that, Eq. (|2.16[ ) takes on the form 

!Lp(n,t) = ^{[w-{n)-w+{nj\ P(n,t)} + ~ {[«T(n) + w+(n)} P(n,t)} (2.25) 

where [w + (n) — u>~(n)] and [w~(n) + u; + (n)] /2 play the roles of drift velocity and diffusion coefficient, respectively. 

To derive the proper boundary conditions we recall that at n = 1 the ME has a different form, given by Eq. (|2.18p . 
It is convenient to require that Eq. (|2.15p is still valid with the additional condition 

[w + (n - l)P(n-l,t)-w-(n)P(n,t)] n=1 = (2.26) 

i.e., the transitions between a fictitious state n = and the state n = 1 are also balanced. 

Similarly, to reconcile the equation at n = N, given by Eq. (|2.19|) . with the general ME, Eq. (|2. 15|) . one should 
impose the condition 

[w- (n + l)P(n + 1, t) - w+(n)P(n, t)} n=N = (2.27) 

which again expresses the balance between an artificial state n = N + 1 and the state n = N. 

In order to gain a deeper insight into the boundary conditions given by Eqs. ([2~26]) and (f2T27|) let us represent 
Eq. (j2~H))) in the form 

4j.P(n,t) = A [w-(ri)P{n,t)-w + {n-l)P[n - l,t)] (2.28) 

This representation looks like a discrete version of the continuity equation, stating that the value in the square brackets 
is the probability current (with a negative sign), i.e., 

J(n) = w + {n - l)P(n - l,t) - w~(n)P(n,t) (2.29) 

A comparison of Eq. (|2.29|) with Eqs. (|2 . 26[) and (|2 . 2T[) allows one to conclude that 

J(n = 1) = and J(n = N + 1) = (2.30) 

i.e., one should impose reflecting boundary conditions on both ends of the interval. 
Within the Fokker-Planck formalism the probability current has the form 

1 d 

J(n) = \w+ (n) - w- (n)l Pin, t) - - — {\w + in) + in)} Pin, t)\ (2.31) 
L J 2 on 

Thus the Fokker-Planck formalism makes it possible to map the strong adsorption case onto a one-dimensional 
random walk problem with drift and diffusion coefficients given in terms of rate constants, Eq. (|2.25p . While such a 
description provides physical insight into the problem, from the viewpoint of numerics it is much easier to deal with 
the ME discrete set Eqs. (j2~T5|) . (|2~18)) and ([2TT9]) . We will discuss the results of this solution in Sec. II D. 



C. Train distribution 



Our MC-simulation results show that the distribution of loops in case of strong physisorption is mainly dominated 
by the shortest loops of length unity. These loops can be considered as defects during the process of zipping. Moreover, 
this distribution sets on much faster than the time for complete adsorption. Thus one may consider the total number 
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of the adsorbed monomers n(t) as a slow variable in comparison to the number of defects (or loops of length unity). 
The adsorbed monomers can be seen as an array of trains, separated by an equilibrium number of defects (see Fig. 
[2}d). The partition function of this one-dimensional array can be determined rigorously (see Appendix A). Thus, one 
derives an expression for the train distribution function 



D(h,t) 



1 



exp 



h av (t) 



(2.32) 



where h av (t) is the average train length. Eq. (12.32[) is nothing but the Flory-Schulz distribution which usually governs 
the molecular weight distribution in equilibrium polymerization of a broad class of systems, referred to as living 
polymers (2~ij]. 



D. Results from the ME numerical solution 



The set of ordinary differential equations (|2.15[) , (|2.18[) and (|2.19p with the initial condition, Eq. (|2.20p . has been 
solved numerically in this investigation. Typically, we use a chain length N = 32, the total time interval takes 300 
units of the the elementary time tq = a Co/hsT, the sticking energy was chosen (in units of fceT) as e = 4.0, whereas 
the entropy loss ln(/i 3 /^ 2 ) = hi2. Figure [3] demonstrates the result of this solution. 




FIG. 3: Adsorbed monomer number distribution function P(n,t) (a) and its isolines as a 2-dimensional log- log plot (b). The 
variation of the distribution maximum, n max (t), is a straight (dashed) line with slope 0.63. 

As it can be seen from FigJ31 the adsorbtion kinetics follows indeed the drift - diffusion picture. The initial 
distribution is very narrow: the adsorption starts with n(0) = I as a grafted chain configuration. As time goes by, the 
distribution maximum moves to larger adsorbed monomer numbers and the distribution itself broadens. Eventually, 
the random process hits the boundary n = N and stays there due to drift and the reflecting boundary conditions. 
As a result, the final distribution is a very narrow again, and is concentrated around the boundary n — N . ft is of 
interest that in the double logarithmic coordinates the distribution maximum follows a straight line (cf. Fig. [3] right 
panel) which reveals a clear scaling law. The first moment n(t) of the distribution function P(n,t) also exhibits well 
expressed scaling behavior, n(t) ~ t ' 66 , as shown in Fig. 2] In the inset we also show the resulting relationship for 
the time of adsorption, t oc A 16 , as expected from Eq. 12.81 Based on the numerical results for P(n,t) and making 
use the relation, Eq. (|2. 10[) . one can calculate the tail distribution function T(l, t) as well. We will discuss this in Sec. 
Ill where we present our MC-results. There it will be seen that our MC-findings are in a good agreement with these 
theoretical predictions. 



III. MONTE CARLO SIMULATION MODEL 



To check the theoretical predictions mentioned in the previous sections we have performed Monte Carlo simulations 
and investigated the adsorption kinetics of a homopolymer, multi-block copolymers, and random copolymers on flat 
surfaces. We have used a coarse grained off-lattice bead spring model[25[ to describe the polymer chains. Our system 
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FIG. 4: The average adsorbed number of monomer vs. time for different chain lengths N. Dashed line denotes the slope, t ' 66 , 
following from Eg. I2.8I In the inset we show the resulting scaling of the adsorption time with chain length, r oc N 1 ' 6 . 

consists of a single chain tethered at one end to a flat structureless surface. There are two kinds of monomers: "A" 
and "B", of which only the "A" type feels an attraction to the surface. The surface interaction of the "A" type 
monomers is described by a square well potential U w (z) = e for z < 6 and U w (z) = otherwise. Here e/fc^T is varied 
from 2.5 to 10.0. The effective bonded interaction is described by the FENE (finitely extensible nonlinear elastic) 
potential. 



Ufene 



-K(l-l ) 2 \n 



1 



l-lo 



lo 



(3.1) 



with K = 20, l m ax = 1, Zq = 0.7, l min = 0.4 

The nonbonded interactions are described by the Morse potential. 

U M (r) 



CM 



exp(-2a(r - r mm )) - 2exp(-a(r - r mm )) 



(3.2) 



with a = 24, r mm = 0.8, e M /k B T = 1. 

We use periodic boundary conditions in the x — y directions and impenetrable walls in the z direction. We 
have studied polymer chains of lengths 32, 64, 128, 256 and 512. Apart from homopolymers, we have also studied 
copolymer chains with block size M = 1 16 and random copolymers (with a fraction of attractive monomers, 
p = 0.25, 0.5, 0.75). The size of the box was 64 x 64 x 64 in all cases except for the 512 chains where we used a larger 
box size of 128 x 128 x 128. The standard Metropolis algorithm was employed to govern the moves with self avoidance 
automatically incorporated in the potentials. In each Monte Carlo update, a monomer was chosen at random and a 
random displacement attempted with Ax, Ay, Az chosen uniformly from the interval —0.5 < Ax, Ay, Az < 0.5. The 
transition probability for the attempted move was calculated from the change AU of the potential energies before and 
after the move as W = exp(— AU /ksT). As for a standard Metropolis algorithm, the attempted move was accepted if 
W exceeds a random number uniformly distributed in the interval [0, 1). A Monte Carlo Step (MCS) is elapsed when 
all N monomers of the chain are selected at random, and given the chance to perform an elementary move. Before 
the surface adsorption potential is switched on, the polymer chain is equilibrated by the MC method for a period of 
about 10 6 MCS (depending on the chain length N this period is varied) whereupon one performs 200 measurement 
runs, each of length 8 x 10 6 MCS. In the case of random copolymers, for a given composition, i.e., percentage p of the 
A— monomers, we create a new polymer chain in the beginning of the simulation run by means of a randomly chosen 
sequence of segments. This chain is then sampled during the course of the run, and replaced by a new sequence in 
the beginning of the next run. 



IV. MONTE CARLO SIMULATION RESULTS 



We present here the main results from the computer simulation of the adsorption kinetics and compare them to 
those from the solution of the Master Equation, Eqs. (|2.15|) . (I2.18P and (|2.19|) . validating thus the theoretical picture 
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of Section HT1 



A. Order Parameter Kinetics - homopolymers 

In Fig. [5^i we show the adsorption time transients which describe the time variation of the order parameter n(t)/N 
(the fraction of adsorbed segments) for homopolymer chains of different length N and strong adhesion e/fceT — 4.0. 
Evidently, in log-log coordinates these transients appear as straight lines, suggesting that the time evolution of the 
adsorption process is governed by a power law. As the chain length N is increased, the slope of the curves grows 
steadily, and for length N — 256 it is equal to » 0.56. This value is close to the theoretically expected slope of 
(1 + v)~ x w 0.62 - cf. Eq. 12. 8[ and for even longer lengths of the polymers would most probably be observed. The 
total time r it takes a polymer chain to be fully adsorbed can be determined from the intersection of the respective 
late time plateau of each transient with the straight line tangent to this transient. Thus one may check the scaling 
of r with polymer length N. In the inset to Fig. [5^i we show the observed scaling of the adsorption time with chain 
length, r oc N a whereby the observed power a rs 1.51 is again somewhat smaller than the expected one 1 + v rj 1.59. 
This small discrepancy is most probably due to finite-size effects too. 

Fig. [5Jd presents the adsorption transients for a chain of constant length, N — 256, for different strength of the 
surface potential. Evidently, as the surface potential gets stronger, the final (equilibrium) values of the transients at 
late times t — * oo grow while the curves are horizontally shifted to shorter times. Notwithstanding, the slope of the 
n(t) curves remains unchanged when e/fegT is varied, suggesting that the kinetics of the process is well described by 
the assumed zipping mechanism. 




FIG. 5: (a) Time evolution of the order parameter (fraction of adsorbed segments) for four different chain lengths TV = 
32, 64, 128, and 256 at surface potential e/fcsT = 4.0. The slope of the N = 256-curve is 0.56. The inset shows the scaling of 
the adsorption time with chain length, r oc N 1 ' 51 . The time r is determined from the intersection point of the late time plateau 
with the tangent t 0,56 to the respective n(t)-curve. (b) Adsorption kinetics for different strengths e of the surface potential. 
The variation of the plateau height (i.e., the fraction of adsorbed monomers at equilibrium) with e is depicted in the upper 

inset where the solid line nt^oo = 1 — 5exp ^— FbtJ describes the equilibrium number of defects (vacancies). The lower inset 
shows a collapse of the adsorption transients on a single 'master curve', if the time axis is rescaled appropriately. 



The changing plateau height may readily be understood as reflecting the correction in the equilibrium fraction of 
adsorbed monomers due to the presence of defects (vacancies) for any given value of e/ksT. This is demonstrated in 
the upper left inset in Fig. 03 where the observed plateau values are shown to be perfectly described by the expression 

nt—>oc — 1 — 5exp (^—jrpf^J under the assumption that the probability of a monomer to desorb from the surface (and 

create a vacancy in the train) is determined by the Boltzmann factor exp -^f- Evidently, the factor of 5 in front of 
the exponent yields the entropic gain in free energy when an adsorbed monomer detaches from the surface while its 
nearest neighbors still stick to it. 

The second inset in Fig. [SJd shows that the adsorption time transients collapse on a master curve, if one rescales the 
time axis appropriately. Note that for a very strong potential, c/KbT = 10.0, the corresponding transient deviates 
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somewhat from the master curve since the establishment of local equilibrium (which we assumed in the theory to 
happen much faster than the adsorption process itself) is hampered. Also the transient for e/kgT = 2.5 (not shown 
in this inset) was found not to fit into the master curve since this strength is close to that of the critical threshold for 
adsorption, the attraction to the surface is comparatively weak and zipping is not the adequate mechanism. For the 
transients which do collapse on a master curve, however, one may view the rescaling of the time axis in Fig. [SJd by 
the expression t — > t[l — 13.7 exp jr^] as a direct confirmation of Eq. 12.231 where the time variable t may be rescaled 
with the driving force of the process (i.e., with the expression in square brackets). The factor w 13.7 gives then the 
ratio 1^3/1^2 of the effective coordination numbers in 3- and 2-dimensions of a polymer chain with excluded volume 
interactions. and \i<i are model-dependent and characterize, therefore, our off-lattice model. 

B. Order Parameter Kinetics - regular and random copolymers 

In Fig. [6] we examine the adsorption kinetics for the case of regular block copolymers with block size M - Fig. [6^,, and 
for random copolymers - Fig.[6jD, bearing in mind that the zipping mechanism, assumed in our theoretical treatment, 
is by no means self-evident when the file of sticking A-monomers is interrupted by neutral B segments. It becomes 
evident from Fig. [5^, however, that, except for a characteristic 'shoulder' in the adsorption transients, the power-law 
character of the order parameter variation with time remains unchanged. Evidently, only the first shoulder in the 
adsorption transient is well expressed while the subsequent ones disappear agaist the background of much larger time 
scales in the log-log representation of Fig. HJi. If, however, one monitors the adsorption of only a single adsorption 
event with time then one observes in normal coordinates a series of such shoulders like a 'staircase' in N a d s (t) (not 
shown here). 

The variation of the power exponent, a, with block length M , where a describes the scaling of the total adsorption 
time with polymer size N , t ox N a , is displayed in the inset right. Evidently, a declines as the block size is increased. 
This finding appears surprising at first sight, since it goes against the general trend of regular multiblock copolymers 
resembling more and more homopolymers (with a = 1 + v for the latter), as the block size M — > 00. Moreover, it 
would imply shorter adsorption times for smaller block size, M — > 1, although the shoulder length visibly grows with 
growing M - see Fig. In fact, however, as one may readily verify from Fig. the transients are systematically 
shifted to longer times (i.e., the total adsorption takes longer) due to a growing prefactor for M — » 1 which does not 
alter the scaling relationship r oc N a . One may thus conclude that the frequent disruption of the zipping process 
for smaller blocks M slows down the overall adsorption process (a transient 'staircase' with numerous short steps) in 
comparison to chains with larger M where the zipping mechanism is fast (a 'staircase' with few longer steps). 




t [MCS] t[MCS] 



FIG. 6: (a) Number of adsorbed segments, N a d s (t)), versus time t for regular AB-copolymers with block size M = 1 -5- 64 and 
length N = 256. For comparison, the transient of a homopolymer is shown by a solid line too. The time interval, taken by 
the initial "shoulder", is shown in the upper left inset. The lower inset displays the variation of the scaling exponent, a, for 
the time of adsorption r oc N a versus block length relationship, (b) The same as in (a) but for random copolymers of length 
N — 256 and different composition p — 0.25, 0.5, 0.75. For p = 1 one has the case of a homopolymer. The inset shows the 
variation of a with p. 
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The characteristic shoulder in the adsorption transients of regular multiblock copolymers manifests itself in the early 
stage of adsorption and lasts progressively longer when M grows. We interpret the temporal length of this shoulder 
with the time it takes for a segment from the second adsorptive A-block in the polymer chain to be eventually captured 
by the attractive surface, once the first A-block has been entirely adsorbed. For sufficiently large blocks one would 
therefore expect that this time interval, t s , associated with the capture event, will scale as the Rouse time, M 1+2ly , 
of a non-adsorbing tethered chain of length M. The observed r s versus M relationship has been shown in the upper 
left inset in Fig. [5^. The slope of w 1.49 is less that the Rouse time scaling exponent, 2.18, which one may attribute 
to the rather small values of the block length M that were accessible in our simulation. One should also allow for 
scatter in the end time of the shoulder due to the mismatch in the capture times of all the successive A-blocks in the 
course of our statistical everaging over many chains during the computer experiment. 

In the case of random copolymers, Fig.[5jD, the observed adsorption transients resemble largely those of a homopoly- 
mer chain with the same number of beads again, apart from the expected difference in the plateau height which is 
determined by the equilibrium number of adsorbed monomers. One should note, however, that a rescaling of the 
vertical axis with the fraction of sticking monomers, p, does not lead to coinciding plateau height - evidently the loops 
whose size also depends on p also affect the equilibrium number of adsorbed monomers. The variation of the observed 
scaling exponent a with composition p is shown in the inset to Fig. [5Jd wherefrom one gets a « 1.6. Note that this 
value is considerably lower than the power of 2.24 which has been observed earlier [hi ], however, for very short chains 
with only 10 sticking beads. One may conclude that even for random copolymer adsorption the typical time of the 
process scales as r cx N a , as observed for homo - and regular block copolymers. It is conceivable, therefore, that an 
effective zipping mechanism in terms of renormalized segments, that is, segments consisting of an A and B diblock 
unit of length 2M for regular multiblock copolymers provides an adequate notion of the way the adsorption kinetics 
may be treated even in such more complicated cases. For random copolymers the role of the block length M would 
then be played by the typical correlation length. 

C. Probability Distribution Functions 

The time evolution in the corresponding Probability Distribution Functions (PDF) of all the trains, loops and 
tails of adsorbed polymers provides a lot of information and insight in the kinetics of the adsorption process. In the 
Appendix we have derived theoretically the expected train distribution under the assumption that local equilibrium 
of loops of unit length is established much faster than the characteristic time of adsorption itself. The resulting 
distribution of possible train lengths is shown to be exponential, in close analogy to that of living polymers [24j]. In 
Fig. [7^ we plot the observed PDF of train lengths for a chain with TV = 256 at two strengths e/ksT of the adsorption 
potential. When scaled with the mean train length h av (t) = (h(t)), at time t, in both cases for e/ksT = 3.0 and 5.0 
one finds an almost perfect straight line in semi-log coordinates, as predicted by Eq. 12.321 



One may thus conclude that the PDF for train lengths preserves its exponential form during the course of the 
adsorption process, validating thus the conjecture of rapid local equilibrium. The latter, however, is somewhat 
violated for the case of rather strong adsorption - e/ksT — 5.0 - shown in Fig. [7^ which is manifested by the increased 
scatter of data at late times when the adsorption process overtakes to some extent the relaxation kinetics on the 
surface. The PDF of loops W(k, t) at different times after the onset of adsorption is shown in Fig. [TJd. Evidently, the 
distribution is sharply peaked at size one whereas less than the remaining 20% of the loops are of size two. Thus the 
loops can be viewed as single thermally activated defects (vacancies) consisting of a desorbed single bead with both 
of its nearest neighbors still attached to the adsorption plane. As the inset in Fig. [7Jd indicates, the PDF of loops is 
also described by an exponential function. The PDFs for loops at different time collapse on a master curve, if scaled 
appropriately with the instantaneous order parameter n(t)/N. 



Eventually, in Fig. [5^ we present the observed PDF of tails for different times t after the start of adsorption, and 
compare the simulation results with those from the numeric solution for T(l, t) according to Eq. 12.101 One may readily 
verify from Fig. [5] that the similarity between simulational and theoretic results is really strong. In both cases one 
starts at t = 1 with a strongly peaked PDF at the full tail length l(t = 1) = N. As time goes by, the distribution 
becomes broader and its maximum shifts to smaller values. At late times the moving peak shrinks again and the 
tail either vanishes, or reduces to a size of single segment which is expressed by the sharp peak at the origin of the 
abscissa. 
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FIG. 7: (a) Distribution of train lengths during the adsorption process of a homopolymer chain with N = 256 at two strengths 
of the adsorption potential e, shown in semi-log coordinates. PDFs for different times (in units of 10 5 MCS) collapse on master 
curves when rescaled by the mean train length h av (t). (b) Distribution of loop lengths W(k,t) for N = 256 and e/ksT = 4.0 
d"™~ „Ac.^„±;„„ t„ +v,„ +v,„ Dr-iT? ;„ normalized by n(t) and shown to be a straight line in log-log 

c 




FIG. 8: (a) Distribution of tail size for different times (in units of 10 5 MCS) during the polymer chain adsorption for a chain 
with N = 256 at e/ksT — 4.0. (b) The same as in (a) as derived from the solution of the ME for chain length N = 32. For 
better visibility the time slices for t = 1, 5, 30 100, 150, 200, and 300 are shifted along the time axis and arranged such that 
the initial distribution for t = 1 is represented by the most distant slice. 



V. CONCLUSIONS 



In this study we examine the adsorption kinetics of a single polymer chain on a flat structureless plane in the strong 
physisorption regime. Adopting the stem-flower model for a chain conformation during adsorption, and assuming the 
segment attachment process to follow a "zipping" mechanism, we develop a scaling theory which describes the time 
evolution of the fraction of adsorbed monomers for polymer chains of arbitrary length N at adsorption strength of 
the surface e/fc^T. 

We derive a Master Equation as well as the corresponding Fokker-Planck equation for the time-dependent PDF 
of the number of adsorbed monomers and for the complementary PDF of tails, and define the appropriate reflecting 
boundary conditions. Inherent in this derivation is the assumed condition of detailed balance which makes it possible 
to relate the elementary steps of adsorption/desorption. From the numeric solution of the equivalent discrete set of 
coupled first-order differential equations we find that the growth of the adsorbed fraction of monomers with time is 
governed by a power law, nit) oc 1 1+ 77 , while the typical time of adsorption r scales with the length of the polymer N 
as t a N a with a = 1 + v. The adsorption transients, found in the Monte Carlo simulation are in good agreement 
with these predictions, if one takes into account the finite-size effects due to the finite length of the studied polymer 
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chains. 

We demonstrate also that the height of the long time plateau in the adsorption transients is determined by the 
equilibrium number of vacancies (defects) in the trains of adsorbed monomers. The transients themselves are found 
to collapse on a single master curve, if time is measured in reduced units which scale with the corresponding driving 
force for adsorption as determined by the surface potential e/ksT. 

A deeper insight into the adsorption kinetics is provided by our detailed study of the relevant probability distribu- 
tions of trains, loops and tails during the adsorption. The predicted exponential expression for the PDF of trains is 
in a very good agreement with our simulational findings. The loops in the strong physisorption regime are observed 
to reduce to occasional desorbed segments (vacancies) which play little role in the dominating picture of trains and 
tails. The PDFs of the latter are found from the simulation data to present a shape which is fully consistent with that 
of the theoretic treatment. It should be noted also that for chemisorption, a monomer adsorption event involves a 
significant local activation barrier [l7])[ll|]. In this so-called "accelerated zipping" regime, the loops formation disrupts 
the adsorption process and the corresponding dynamics differs significantly from the one investigated in this paper. 

Eventually, in the case of regular multiblock and random copolymers we find that the adsorption kinetics strongly 
resembles that of homopolymers. The observed deviations from the latter suggest plausible interpretations in terms 
of polymer dynamics, however, it is clear that additional investigations will be warranted before a complete picture 
of the adsorption kinetics in this case is established too. 
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APPENDIX A: DERIVATION OF TRAIN DISTRIBUTION 



The partition function of an one-dimensional array of p + 1 trains, separated by p defects, has the following form 

$[n(t),p] = J ...J dx\...dx v 

Q<x 1 <X2...x p <n(t) 
n(t) n(t) n(t) 

dx x J dx 2 ... J dx p = [n(t)] p (Al) 

xi i p -i 

where n{t) is the total number of adsorbed monomers at time t. 

Consider now the the distribution of an arbitrary train h s+ i = x s+ i — x s . In order to find it, one should carry out 
the integration in Eq. (|A1[) over all ^-coordinates except x s and x s +\. In result of the integration one gets 

$xsx.+i[n(t),p] dx s dx s+1 = — - ~ _ ^| x t^ 1 K f ) _ ^s+i] P ~ S ~ 1 dx s dx s+1 (A2) 

where Eq. ljAip has been used separately for the intervals [0, x s ] and [x s+ i,n(t)]. 

The distribution of the train length, h s+ i = x s+ i — x s , follows immediately from Eq. (|A2[) after integrating over x s , 
i.e. 

n(i)—h s+ i 

*h. +l [n(t),p] - (s _ 1)! J_ s _ 1) , J X s ' 1 [n(t) - h s+1 - x^- 1 dx s (A3) 

o 

By the substitution, y = x s /[n(t) — h s+ i], in the integral of ea. (|A3p one arrives at the result 



l 



**^[»(t),ri = , [n( *L*' +ir 1 \, Iv-Hi-vY"- 1 *, 



(s-l)l(p-s-l)! 





(A4) 
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where one has used Jq y s ~ 1 (l — y) p ~ s ~ 1 dy = (s — l)!(p — s — l)!/(p — 1)!- The result in Eq. (IA4j) does not depend on 
the consecutive number of the train, as expected. The normalized probability to find a train of the length h at time 
t is given by 



D(h,t) = 



*h[n(t),p] = pi [n(t) - h] 
t>[n{t),p] (p-1)! [n(t)f 



p-i 



P 
n(t) 



h 



p-1 



n{t) 



exp 



n(t) 



(A5) 



where one uses Eqs. dAll) and (I A4|) as well as the conditions p ^> 1 and h/n(t) -C 1. Taking into account that the 
average train length h av (t) = n(t)/p, the last expression results in Eq. (|2.32p . 
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